Dual domestications and origin of traits in grapevine evolution.

We elucidate grapevine evolution and domestication histories with 3525 cultivated and wild accessions worldwide. In the Pleistocene, harsh climate drove the separation of wild grape ecotypes caused by continuous habitat fragmentation. Then, domestication occurred concurrently about 11,000 years ago in Western Asia and the Caucasus to yield table and wine grapevines. The Western Asia domesticates dispersed into Europe with early farmers, introgressed with ancient wild western ecotypes, and subsequently diversified along human migration trails into muscat and unique western wine grape ancestries by the late Neolithic. Analyses of domestication traits also reveal new insights into selection for berry palatability, hermaphroditism, muscat flavor, and berry skin color. These data demonstrate the role of the grapevines in the early inception of agriculture across Eurasia.

Interest of phenomic prediction as an alternative to genomic prediction in grapevine.

BACKGROUND: Phenomic prediction has been defined as an alternative to genomic prediction by using spectra instead of molecular markers. A reflectance spectrum provides information on the biochemical composition within a tissue, itself being under genetic determinism. Thus, a relationship matrix built from spectra could potentially capture genetic signal. This new methodology has been mainly applied in several annual crop species but little is known so far about its interest in perennial species. Besides, phenomic prediction has only been tested for a restricted set of traits, mainly related to yield or phenology. This study aims at applying phenomic prediction for the first time in grapevine, using spectra collected on two tissues and over two consecutive years, on two populations and for 15 traits, related to berry composition, phenology, morphological and vigour. A major novelty of this study was to collect spectra and phenotypes several years apart from each other. First, we characterized the genetic signal in spectra and under which condition it could be maximized, then phenomic predictive ability was compared to genomic predictive ability. RESULTS: For the first time, we showed that the similarity between spectra and genomic relationship matrices was stable across tissues or years, but variable across populations, with co-inertia around 0.3 and 0.6 for diversity panel and half-diallel populations, respectively. Applying a mixed model on spectra data increased phenomic predictive ability, while using spectra collected on wood or leaves from one year or another had less impact. Differences between populations were also observed for predictive ability of phenomic prediction, with an average of 0.27 for the diversity panel and 0.35 for the half-diallel. For both populations, a significant positive correlation was found across traits between predictive ability of genomic and phenomic predictions. CONCLUSION: NIRS is a new low-cost alternative to genotyping for predicting complex traits in perennial species such as grapevine. Having spectra and phenotypes from different years allowed us to exclude genotype-by-environment interactions and confirms that phenomic prediction can rely only on genetics.

A genome-wide association and prediction study in grapevine deciphers the genetic architecture of multiple traits and identifies genes under many new QTLs.

To cope with the challenges facing agriculture, speeding-up breeding programs is a worthy endeavor, especially for perennial species such as grapevine, but requires understanding the genetic architecture of target traits. To go beyond the mapping of quantitative trait loci in bi-parental crosses, we exploited a diversity panel of 279 Vitis vinifera L. cultivars planted in 5 blocks in the vineyard. This panel was phenotyped over several years for 127 traits including yield components, organic acids, aroma precursors, polyphenols, and a water stress indicator. The panel was genotyped for 63k single nucleotide polymorphisms by combining an 18K microarray and genotyping-by-sequencing. The experimental design allowed to reliably assess the genotypic values for most traits. Marker densification via genotyping-by-sequencing markedly increased the proportion of genetic variance explained by single nucleotide polymorphisms, and 2 multi-single nucleotide polymorphism models identified quantitative trait loci not found by a single nucleotide polymorphism-by-single nucleotide polymorphism model. Overall, 489 reliable quantitative trait loci were detected for 41% more response variables than by a single nucleotide polymorphism-by-single nucleotide polymorphism model with microarray-only single nucleotide polymorphisms, many new ones compared with the results from bi-parental crosses. A prediction accuracy higher than 0.42 was obtained for 50% of the response variables. Our overall approach as well as quantitative trait locus and prediction results provide insights into the genetic architecture of target traits. New candidate genes and the application into breeding are discussed.

Across-population genomic prediction in grapevine opens up promising prospects for breeding.

Crop breeding involves two selection steps: choosing progenitors and selecting individuals within progenies. Genomic prediction, based on genome-wide marker estimation of genetic values, could facilitate these steps. However, its potential usefulness in grapevine (Vitis vinifera L.) has only been evaluated in non-breeding contexts mainly through cross-validation within a single population. We tested across-population genomic prediction in a more realistic breeding configuration, from a diversity panel to ten bi-parental crosses connected within a half-diallel mating design. Prediction quality was evaluated over 15 traits of interest (related to yield, berry composition, phenology and vigour), for both the average genetic value of each cross (cross mean) and the genetic values of individuals within each cross (individual values). Genomic prediction in these conditions was found useful: for cross mean, average per-trait predictive ability was 0.6, while per-cross predictive ability was halved on average, but reached a maximum of 0.7. Mean predictive ability for individual values within crosses was 0.26, about half the within-half-diallel value taken as a reference. For some traits and/or crosses, these across-population predictive ability values are promising for implementing genomic selection in grapevine breeding. This study also provided key insights on variables affecting predictive ability. Per-cross predictive ability was well predicted by genetic distance between parents and when this predictive ability was below 0.6, it was improved by training set optimization. For individual values, predictive ability mostly depended on trait-related variables (magnitude of the cross effect and heritability). These results will greatly help designing grapevine breeding programs assisted by genomic prediction.

Harnessing multivariate, penalized regression methods for genomic prediction and QTL detection of drought-related traits in grapevine.

Viticulture has to cope with climate change and to decrease pesticide inputs, while maintaining yield and wine quality. Breeding is a key lever to meet this challenge, and genomic prediction a promising tool to accelerate breeding programs. Multivariate methods are potentially more accurate than univariate ones. Moreover, some prediction methods also provide marker selection, thus allowing quantitative trait loci (QTLs) detection and the identification of positional candidate genes. To study both genomic prediction and QTL detection for drought-related traits in grapevine, we applied several methods, interval mapping (IM) as well as univariate and multivariate penalized regression, in a bi-parental progeny. With a dense genetic map, we simulated two traits under four QTL configurations. The penalized regression method Elastic Net (EN) for genomic prediction, and controlling the marginal False Discovery Rate on EN selected markers to prioritize the QTLs. Indeed, penalized methods were more powerful than IM for QTL detection across various genetic architectures. Multivariate prediction did not perform better than its univariate counterpart, despite strong genetic correlation between traits. Using 14 traits measured in semi-controlled conditions under different watering conditions, penalized regression methods proved very efficient for intra-population prediction whatever the genetic architecture of the trait, with predictive abilities reaching 0.68. Compared to a previous study on the same traits, these methods applied on a denser map found new QTLs controlling traits linked to drought tolerance and provided relevant candidate genes. Overall, these findings provide a strong evidence base for implementing genomic prediction in grapevine breeding.

Genetic diversity and population structure in Vitis species illustrate phylogeographic patterns in eastern North America.

Geographical distribution and diversity of current plant species have been strongly shaped by climatic oscillations during the Quaternary. Analysing the resulting divergence among species and differentiation within species is crucial to understand the evolution of taxa like the Vitis genus, which provides very useful genetic resources for grapevine improvement and might reveal original recolonization patterns due to growth habit and dispersal mode. Here, we studied the genetic structure in natural populations of three species from eastern North America: Vitis aestivalis, V. cinerea and V. riparia using different marker types. Vitis aestivalis and V. cinerea showed higher diversity than V. riparia. The two former species are less differentiated, confirming an earlier divergence of V. riparia. V. aestivalis and V. riparia exhibited different genetic groups on both sides of the Appalachian Mountains that could mirror different recolonization routes from southern refugia. Genetic structure was stronger in V. cinerea, for which two varieties (var. berlandieri and var. cinerea) are morphologically recognized. Our results confirm this distinction and suggest the existence of three other lineages within var. cinerea. These discontinuities appear linked to adaptation of var. berlandieri to dry and limy areas of Texas and partially to the Mississippi River Valley. Rapid range expansions from refugia upon climate warming are also suggested by the low linkage disequilibrium values observed. Furthermore, large variation for downy mildew resistance was observed in the three species. Our findings appeared consistent with the vegetation history of eastern North America.

Beyond Seek and Destroy: how to Generate Allelic Series Using Genome Editing Tools.

Genome editing tools have greatly facilitated the functional analysis of genes of interest by targeted mutagenesis. Many usable genome editing tools, including different site-specific nucleases and editor databases that allow single-nucleotide polymorphisms (SNPs) to be introduced at a given site, are now available. These tools can be used to generate high allelic diversity at a given locus to facilitate gene function studies, including examining the role of a specific protein domain or a single amino acid. We compared the effects, efficiencies and mutation types generated by our LbCPF1, SpCAS9 and base editor (BECAS9) constructs for the OsCAO1 gene. SpCAS9 and LbCPF1 have similar efficiencies in generating mutations but differ in the types of mutations induced, with the majority of changes being single-nucleotide insertions and short deletions for SpCAS9 and LbCPF1, respectively. The proportions of heterozygotes also differed, representing a majority in our LbCPF1, while with SpCAS9, we obtained a large number of biallelic mutants. Finally, we demonstrated that it is possible to specifically introduce stop codons using the BECAS9 with an acceptable efficiency of approximately 20%. Based on these results, a rational choice among these three alternatives may be made depending on the type of mutation that one wishes to introduce, the three systems being complementary. SpCAS9 remains the best choice to generate KO mutations in primary transformants, while if the desired gene mutation interferes with regeneration or viability, the use of our LbCPF1 construction will be preferred, because it produces mainly heterozygotes. LbCPF1 has been described in other studies as being as effective as SpCAS9 in generating homozygous and biallelic mutations. It will remain to be clarified in the future, whether the different LbCFP1 constructions have different efficiencies and determine the origin of these differences. Finally, if one wishes to specifically introduce stop codons, BECAS9 is a viable and efficient alternative, although it has a lower efficiency than SpCAS9 and LbCPF1 for creating KO mutations.

Understanding the phyllosphere microbiome assemblage in grape species (Vitaceae) with amplicon sequence data structures.

Impacts of plant genotype on microbial assemblage in the phyllosphere (above-ground parts of plants, which predominantly consists of the set of photosynthetic leaves) of Vitis vinifera cultivars have been studied previously but the impact of grape species (under the grape family Vitaceae) was never investigated. Considering the fact, that the phyllosphere microbiome may have profound effects on host plant health and its performance traits, studying the impact of grape species in microbial taxa structuring in the phyllosphere could be of crucial importance. We performed 16S and ITS profiling (for bacteria and fungi respectively) to access genus level characterization of the microflora present in the leaf phyllosphere of five species within this plant family, sampled in two successive years from the repository situated in the Mediterranean. We also performed α and ß-diversity analyses with robust statistical estimates to test the impacts of grape species and growing year, over a two-year period. Our results indicated the presence of complex microbial diversity and assemblages in the phyllosphere with a significant effect of both factors (grape species and growing year), the latter effect is being more pronounced. We also compared separate normalization methods for high-throughput microbiome data-sets followed by differential taxa abundance analyses. The results suggested the predominance of a particular normalization method over others. This also indicated the need for more robust normalization methods to study the differential taxa abundance among groups in microbiome research.

Palaeogenomic insights into the origins of French grapevine diversity.

The Eurasian grapevine (Vitis vinifera) has long been important for wine production as well as being a food source. Despite being clonally propagated, modern cultivars exhibit great morphological and genetic diversity, with thousands of varieties described in historic and contemporaneous records. Through historical accounts, some varieties can be traced to the Middle Ages, but the genetic relationships between ancient and modern vines remain unknown. We present target-enriched genome-wide sequencing data from 28 archaeological grape seeds dating to the Iron Age, Roman era and medieval period. When compared with domesticated and wild accessions, we found that the archaeological samples were closely related to western European cultivars used for winemaking today. We identified seeds with identical genetic signatures present at different Roman sites, as well as seeds sharing parent-offspring relationships with varieties grown today. Furthermore, we discovered that one seed dated to ~1100 CE was a genetic match to 'Savagnin Blanc', providing evidence for 900 years of uninterrupted vegetative propagation.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population.

KEY MESSAGE: A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

Genotype-Environment Interaction Shapes the Microbial Assemblage in Grapevine's Phyllosphere and Carposphere: An NGS Approach.

Plant surface or phyllosphere is the habitat of hyperdiverse microbial communities and it is always exposed to the fluctuating environmental factors, which is thought to be one of the potential drivers of microbial community structuring. Impact of grapevine genotypes in variable environmental factors (i.e., at different geographic locations) on the phyllosphere has never been studied and is the main objective of this report. Using high throughput short amplicon sequencing of 16S rRNA genes and internal transcribed spacer (ITS), we analyzed the impacts of genotypes of Vitis Vinifera (coming from three genetic pool), on the microbial (bacterial and fungal) assemblage in the phyllosphere. First, we performed the analysis of the phyllosphere microbiome while using fifteen genotypes that were chosen to maximize intra-specific diversity and grown in two Mediterranean vineyards. Then, the same analysis was performed on five commercially important varieties of Vitis vinifera that were sampled from three different French agro-climatic zones (or terroir: a combination of climate, soils, and human practices). Our study revealed that, at a particular geographic location, genotypes have an impact on microbial assemblage in the phyllosphere and carposphere of leaf and fruit (or berries), respectively, which is more prominent on the carposphere but the effect of terroir was much stronger than the genotype when the leaf phyllosphere of five grapevine varieties grown in different agro-climatic zones was compared. Impacts of the season and exterior plant organs (leaf and berries) on microbial taxa structuring in the phyllosphere was also assessed and presented in this report.

Early Neolithic wine of Georgia in the South Caucasus.

Chemical analyses of ancient organic compounds absorbed into the pottery fabrics from sites in Georgia in the South Caucasus region, dating to the early Neolithic period (ca. 6,000-5,000 BC), provide the earliest biomolecular archaeological evidence for grape wine and viniculture from the Near East, at ca. 6,000-5,800 BC. The chemical findings are corroborated by climatic and environmental reconstruction, together with archaeobotanical evidence, including grape pollen, starch, and epidermal remains associated with a jar of similar type and date. The very large-capacity jars, some of the earliest pottery made in the Near East, probably served as combination fermentation, aging, and serving vessels. They are the most numerous pottery type at many sites comprising the so-called "Shulaveri-Shomutepe Culture" of the Neolithic period, which extends into western Azerbaijan and northern Armenia. The discovery of early sixth millennium BC grape wine in this region is crucial to the later history of wine in Europe and the rest of the world.

Genetic diversity, linkage disequilibrium and power of a large grapevine (Vitis vinifera L) diversity panel newly designed for association studies.

BACKGROUND: As for many crops, new high-quality grapevine varieties requiring less pesticide and adapted to climate change are needed. In perennial species, breeding is a long process which can be speeded up by gaining knowledge about quantitative trait loci linked to agronomic traits variation. However, due to the long juvenile period of these species, establishing numerous highly recombinant populations for high resolution mapping is both costly and time-consuming. Genome wide association studies in germplasm panels is an alternative method of choice, since it allows identifying the main quantitative trait loci with high resolution by exploiting past recombination events between cultivars. Such studies require adequate panel design to represent most of the available genetic and phenotypic diversity. Assessing linkage disequilibrium extent and panel power is also needed to determine the marker density required for association studies. RESULTS: Starting from the largest grapevine collection worldwide maintained in Vassal (France), we designed a diversity panel of 279 cultivars with limited relatedness, reflecting the low structuration in three genetic pools resulting from different uses (table vs wine) and geographical origin (East vs West), and including the major founders of modern cultivars. With 20 simple sequence repeat markers and five quantitative traits, we showed that our panel adequately captured most of the genetic and phenotypic diversity existing within the entire Vassal collection. To assess linkage disequilibrium extent and panel power, we genotyped single nucleotide polymorphisms: 372 over four genomic regions and 129 distributed over the whole genome. Linkage disequilibrium, measured by correlation corrected for kinship, reached 0.2 for a physical distance between 9 and 458 Kb depending on genetic pool and genomic region, with varying size of linkage disequilibrium blocks. This panel achieved reasonable power to detect associations between traits with high broad-sense heritability (> 0.7) and causal loci with intermediate allelic frequency and strong effect (explaining > 10 % of total variance). CONCLUSIONS: Our association panel constitutes a new, highly valuable resource for genetic association studies in grapevine, and deserves dissemination to diverse field and greenhouse trials to gain more insight into the genetic control of many agronomic traits and their interaction with the environment.

Identification of stable QTLs for vegetative and reproductive traits in the microvine (Vitis vinifera L.) using the 18 K Infinium chip.

BACKGROUND: The increasing temperature associated with climate change impacts grapevine phenology and development with critical effects on grape yield and composition. Plant breeding has the potential to deliver new cultivars with stable yield and quality under warmer climate conditions, but this requires the identification of stable genetic determinants. This study tested the potentialities of the microvine to boost genetics in grapevine. A mapping population of 129 microvines derived from Picovine x Ugni Blanc flb, was genotyped with the Illumina® 18 K SNP (Single Nucleotide Polymorphism) chip. Forty-three vegetative and reproductive traits were phenotyped outdoors over four cropping cycles, and a subset of 22 traits over two cropping cycles in growth rooms with two contrasted temperatures, in order to map stable QTLs (Quantitative Trait Loci). RESULTS: Ten stable QTLs for berry development and quality or leaf area were identified on the parental maps. A new major QTL explaining up to 44 % of total variance of berry weight was identified on chromosome 7 in Ugni Blanc flb, and co-localized with QTLs for seed number (up to 76 % total variance), major berry acids at green lag phase (up to 35 %), and other yield components (up to 25 %). In addition, a minor QTL for leaf area was found on chromosome 4 of the same parent. In contrast, only minor QTLs for berry acidity and leaf area could be found as moderately stable in Picovine. None of the transporters recently identified as mutated in low acidity apples or Cucurbits were included in the several hundreds of candidate genes underlying the above berry QTLs, which could be reduced to a few dozen candidate genes when a priori pertinent biological functions and organ specific expression were considered. CONCLUSIONS: This study combining the use of microvine and a high throughput genotyping technology was innovative for grapevine genetics. It allowed the identification of 10 stable QTLs, including the first berry acidity QTLs reported so far in a Vitis vinifera intra-specific cross. Robustness of a set of QTLs was assessed with respect to temperature variation.

MybA1 gene diversity across the Vitis genus.

The MybA1 gene in the genus Vitis encodes a transcription factor, belonging to the R2R3 Myb family, that controls the last steps in the anthocyanins biosynthesis pathway. Polymorphism within MybA1 has been associated with color variation in berries of V. vinifera and other Vitis species. In this work, we analyzed the sequence variation in MybA1 both in the subg. Muscadinia and in an extended set of Asian, American and European genotypes of subg. Vitis. Our aims were to infer the evolution of this gene during the speciation process and to identify polymorphisms that could potentially generate changes in gene regulation. The results show that MybA1 experienced many insertions and deletions in non-coding regions but also in the third exon sequence. Owing to the larger set of Vitis species compared here, new indels were identified and the origin of previously described indels was reconsidered. A large number of single nucleotide polymorphisms were found in non-coding regions but also in the sequence coding for the R2R3 domain and the C terminal part of the protein. Some of these changes led to amino acid substitutions and therefore could have modified MybA1 protein activity. Bayesian phylogenetic analysis of all polymorphisms did not provide a consensus tree depicting the geographical partitioning of the species but allowed highlighting several species relationships within subgenus Vitis. Finally, the evolutionary events described could be useful to gain more insight into the role of MybA1 for anthocyanin biosynthesis in grapevine.

Genetic diversity of wild and cultivated grapevine accessions from southeast Turkey.

Wild grapevine genetic diversity in southeast Turkey has not been documented to date. In the present work, in order to clarify the relationships between wild and cultivated grape accessions from southeastern Turkey, 22 nuclear and three chloroplast microsatellite loci were used on 21 wild grapevine Vitis vinifera L. ssp. sylvestris (Gmelin) and 13 cultivated grapevine Vitis vinifera ssp. sativa accessions. The number of alleles per SSR locus ranged from 4 (VVIn16) to 20 (VVIv67) and the mean allele number per locus was 10.09. Expected locus heterozygosity ranged from 0.586 (locus VVIb01) to 0.898 (locus (VVIv67)). The three cpSSR molecular markers presented variation in size both in cultivars and in wild Turkish accessions. Two size variants were detected for cpSSR3 (106 and 107 bp) for cpSSR5 (104 and 105 bp), and for cpSSR10 (115 and 116 bp). The six alleles in wild grapevines fell into three haplotypes B, C and D. A genetic structure according to accessions taxonomic status (wild or cultivated) was revealed by UPGMA analysis. This highlighted a clear separation between domesticated and wild accessions in Turkish germplasm. The results pointed out the need to further collect and characterize this wild and cultivated grapevine germplasm.

Genetic variation in a grapevine progeny (Vitis vinifera L. cvs Grenache×Syrah) reveals inconsistencies between maintenance of daytime leaf water potential and response of transpiration rate under drought.

In the face of water stress, plants evolved with different abilities to limit the decrease in leaf water potential, notably in the daytime (ΨM). So-called isohydric species efficiently maintain high ΨM, whereas anisohydric species cannot prevent ΨM from dropping as soil water deficit develops. The genetic and physiological origins of these differences in (an)isohydric behaviours remain to be clarified. This is of particular interest within species such as Vitis vinifera L. where continuous variation in the level of isohydry has been observed among cultivars. With this objective, a 2 year experiment was conducted on the pseudo-F1 progeny from a cross between the two widespread cultivars Syrah and Grenache using a phenotyping platform coupled to a controlled-environment chamber. Potted plants of all the progeny were analysed for ΨM, transpiration rate, and soil-to-leaf hydraulic conductance, under both well-watered and water deficit conditions. A high genetic variability was found for all the above traits. Four quantitative trait loci (QTLs) were detected for ΨM under water deficit conditions, and 28 other QTLs were detected for the different traits in either condition. Genetic variation in ΨM maintenance under water deficit weakly correlated with drought-induced reduction in transpiration rate in the progeny, and QTLs for both traits did not completely co-localize. This indicates that genetic variation in the control of ΨM under water deficit was not due simply to variation in transpiration sensitivity to soil drying. Possible origins of the diversity in (an)isohydric behaviours in grapevine are discussed on the basis of concurrent variations in soil-to-leaf hydraulic conductance and stomatal control of transpiration.

Genome-wide prediction methods in highly diverse and heterozygous species: proof-of-concept through simulation in grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines.

BACKGROUND: In Vitis vinifera L., domestication induced a dramatic change in flower morphology: the wild sylvestris subspecies is dioecious while hermaphroditism is largely predominant in the domesticated subsp. V. v. vinifera. The characterisation of polymorphisms in genes underlying the sex-determining chromosomal region may help clarify the history of domestication in grapevine and the evolution of sex chromosomes in plants. In the genus Vitis, sex determination is putatively controlled by one major locus with three alleles, male M, hermaphrodite H and female F, with an allelic dominance M > H > F. Previous genetic studies located the sex locus on chromosome 2. We used DNA polymorphisms of geographically diverse V. vinifera genotypes to confirm the position of this locus, to characterise the genetic diversity and traces of selection in candidate genes, and to explore the origin of hermaphroditism. RESULTS: In V. v. sylvestris, a sex-determining region of 154.8 kb, also present in other Vitis species, spans less than 1% of chromosome 2. It displays haplotype diversity, linkage disequilibrium and differentiation that typically correspond to a small XY sex-determining region with XY males and XX females. In male alleles, traces of purifying selection were found for a trehalose phosphatase, an exostosin and a WRKY transcription factor, with strikingly low polymorphism levels between distant geographic regions. Both diversity and network analysis revealed that H alleles are more closely related to M than to F alleles. CONCLUSIONS: Hermaphrodite alleles appear to derive from male alleles of wild grapevines, with successive recombination events allowing import of diversity from the X into the Y chromosomal region and slowing down the expansion of the region into a full heteromorphic chromosome. Our data are consistent with multiple domestication events and show traces of introgression from other Asian Vitis species into the cultivated grapevine gene pool.